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^-i' Abstract. Based on a linear and non-linear study of radial pulsations in 

jy^ . the envelopes of classical novae (Schenker 1999), I discuss the results both 

from the point of view of pulsation theory as well as their consequences 
for current nova models. 
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Starting from initially static envelope structures at various stages 
during the decline of a nova outburst, strong 'running wave' instabili- 
ties have been found that rapidly grow into shocks. Improved analytical 
concepts give a new direction to the interpretation of such highly non- 
\ adiabatic radial pulsations. 

Q\ ■ For direct observational confirmation a search for short period vari- 

C> ability in the UV and soft X-ray is suggested during the very late decline 

phase. Speculative consequences for mass loss scenarios in novae due to 
these instabilities will need some more work in the future. 
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Introduction 



Classical nova outbursts take place in mass transferring binaries, when hydrogen- 
rich material is accreted onto a white dwarf (WD). Once enough matter has ac- 
cumulated, nuclear burning sets in explosively causing a complex and dynamical 
initial phase of the nova outburst. Around visual maximum however, a much 
more steady phase begins, where an extremely large ratio of luminosity to mass 
occurs and the structure of the star resembles that of a giant: a nuclear burning 
shell sitting on top of the WD, covered by an extended envelope of very little 
mass (typically < 1O~ 4 M ). 

Schenker (1999) analyses the pulsational stability of such a structure assum- 
ing a static envelope, or (given the ongoing evolution of the nova on time-scales 
of maybe weeks or months) a sequence of quasi-static envelopes. 

2. Models & Results 

The background models for the linear analysis have been constructed by a linear 
series of core plus envelope models, obtained via a shooting method. Figure la 



1 e-mail: kjs@star.le.ac.uk 



1 



2 



Schenker 



a) 



SpT 

19 FO K7 M8 



b) 




15 

.5 20 

B 

\ 

R 30 

50 

100 
500 - 



6.0 5.5 5.0 4.5 4.0 3.5 

log ( T.„ / K ) 





□ 




O* 




OK 


O* 




OK 


Cw 




OK 


0* 




o* 


Ok 




OK 


Ok 






Ok 




OK 


Ok 


□ 


o< 


Ok 




OK 


Ok 




O 


El 






o 


□ 


o* 






o 


o< 




o 


O 




Of 




□ 

. . T , . 




-1.0 -0.5 


0.0 


0.5 1.0 



Figure 1. a) The quasi-evolutionary sequences and the specific 20 
models analysed — b) Complex Eigenspectrum (ctr, cti) of model B 
in different numerical treatment: squares adiabatic, diamonds NAR, 
asterisks non- adiabatic. 



(left panel) shows their position in the HR diagram. Along these curves a radial, 
non-adiabatic linear stability analysis has been performed utilising the Riccati 
method (Gautschy & Glatzel 1990). For the 20 models specifically marked I 
have also studied the non-linear evolution with adaptive-grid radiation hydro- 
code (RHD, Dorfi & Feuchtinger 1991). 

One major result is plotted in Figure lb (right panel), showing the complex 
Eigenspectrum of model B (M = 0.5 Mq, log(T e g/K) = 4.8). Frequencies a are 
plotted normalised with the free-fall time-scale in full non-adiabatic treatment 
and two different approximations. Most striking are the enormous growth rates 
found, of the same order as the free-fall time-scale. This is the case not only in 
model B, but for an extended range of T c gf . In contrast to classical pulsators on 
the instability strip (e.g. RR Lyr stars) here the adiabatic approximation cannot 
reproduce the frequencies at all. However, a highly suitable approximation for 
such a situation is the Non- Adiabatic Reversible approximation (NAR), also 
used by Glatzel (1994) and Saio, Baker, & Gautschy (1998) in their investigations 
into the nature of strange modes. One typical characteristic is the occurrence 
of complex conjugate pairs of Eigenfrequencies as apparent in Figure lb. 

NAR represents the Eigenmodes almost perfectly around the centre of the 
instability region (model B), not only in frequency (Fig. lb) but also in terms of 
the Eigenfunctions. These Eigenfunctions (cf. Fig. 2) represent complex stand- 
ing, but in their projection running waves (in x, I, and d): outwards for the 
unstable, and inwards for the stable mode of each pair. In Figure 2, fully non- 
adiabatic Eigenfunctions are shown by amplitude and phase, so perfectly com- 
plex conjugate pairs of waves running inwards & outwards would have identical 
arguments but negative phase in each of the Eigenfunctions. This is almost 
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Figure 2. Set of fully non-adiabatic, complex Eigenfunctions for the 
unstable mode of model B with (<tr,cji) = (2.92,-0.57), i.e. from the 
5th pair in Fig. lb. In the first 5 panels, the relative Lagrangian per- 
turbations of radial displacement x, luminosity /, density d, pressure 
p, and temperature t are shown split into amplitude and phase angle. 
The last panel shows the phase difference between density and pressure 
perturbation together with the differential work integral. 



fulfilled and a mathematical requirement for NAR solutions, with the exception 
of the luminosity perturbation / which is effectively set to zero in this approxi- 
mation. 



3. Interpretation & Observational Consequences 

An improved local analysis with gradients on such models reveals their nature 
as radiation-modified acoustic modes. It also allows to understand the spectrum 
found (e.g. in Fig. 1) to some degree. 

The complementary analysis using the RHD code provides further charac- 
teristics of the modes and consequences in the non-linear regime. During the 
linear growth, and more so lateron, an additional energy transport is observed, 
as packets of radiation energy are spiralling outwards in density perturbations 
(similar to a "screw conveyor"). This is likely to be responsible for the new 
"equilibrium" structure of background model, which actually forms a curve on 
the HR diagram which can significantly deviate from the initial position. 

As a possibly related phenomena with non-radial behaviour of envelopes 
under strong radiation fields, Shaviv (2001) has studied the application of Super- 
Eddington winds on novae. In his model spatial variations of opacity allow 
a larger radiation flux to pass, whereas in the 1-dimensional case here only 
temporal modulations provide a similar effect. 



4 



Schenker 



Although the project initially attempted to explain oscillation during the 
so-called "transition phase" (Bianchini, Friedjung, & Brinkmann 1992), a more 
probably scenario for direct observational evidence lies in the predicted X-ray 
oscillations shortly before the turn-off of many novae (including the faster ones 
on more massive WDs, cf. Fig. la). 

4. Conclusions 

It has become evident that static equilibrium models of novae after visual max- 
imum are subject to violent pulsational instabilities. This finding (i) may sug- 
gests a possible link between mass loss in novae and pulsations, but (ii) at least 
excludes a purely quasi-static evolution even for the lowest mass novae (possibly 
systems like V723 Cas and HR Del), where common mass loss mechanisms tend 
to fail. For more massive WDs, only the very late phases of the outburst could 
be modelled. While they do also show instabilities, they still might pulsate even 
during the non-static earlier part. 

Classical novae form an excellent laboratory for pulsations under radia- 
tion dominated conditions, so both fields will benefit from their further detailed 
study. 
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Discussion 

A. Cox : Is there any convection included in your calculations? 

K. Schenker : I do not believe that convective effects are very important in this 
case. These envelopes are strongly dominated by radiation, and the radiative flux 
in the static models (which do include convection via MLT) is indeed much higher 
that its convective counterpart. So the linear analysis neglects the contribution 
of the convective flux altogether, which is consistent with the RHD models that 
do not contain convection at all. 

A. Cox : What is the temperature of the pulsating driving region? 

K. Schenker : The driving temperature lies in the range of the 'iron peak' present 
in the new OPAL opacities, around 10 5 2 K. 



